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Abstract 



I Regulatory T cells (Treg) have recently been identified as playing a central role in allergy and during allergen-specific 
immunotherapy. We have extended our previous mathematical model describing the nonlinear dynamics of Thl-Th2 
CN| .regulation by including Treg cells and their major cytokines. We hypothesize that immunotherapy mainly acts on the 
!— i T cell level and that the decisive process can be regarded as a dynamical phenomenon. The model consists of nonlinear 
. ^ idifferential equations which describe the proliferation and mutual suppression of different T cell subsets. The old version 
J>; of the model was based upon the Thl-Th2 paradigm and is successful in describing the "Thl-Th2 switch" which was 
^| considered the decisive event during specific immunotherapy. In recent years, however, the Thl-Th2 paradigm has been 
questioned and therefore, we have investigated a modified model in order to account for the influence of a regulatory T 
r- — i .cell type. We examined the extended model by means of numerical simulations and analytical methods. As the modified 
model is more complex, we had to develop new methods to portray its characteristics. The concept of stable manifolds 
of fixed points of a stroboscobic map turned out to be especially important. We found that when including regulatory 
q T cells, our model can describe the events in allergen-specific immunotherapy more accurately. Our results suggest that 
• the decisive effect of immunotherapy, the increased proliferation of Treg and suppression of Th2 cells, crucially depends 
F ^ on the administration of high dose injections right before the maintenance phase sets in. Empirical protocols could 
O^.therefore be improved by optimizing this step of therapy. 
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1. Introduction 

T helper cells play a significant role in immune re- 
sponses to allergic substances. There are several subtypes 
of T helper cells which differ in function according to their 
cytokine profiles. Specific Th2 cells are mainly responsible 
for allergic reactions as they can activate the production of 
IgE antibodies by means of which the well known allergic 
symptoms are provoked. The "Thl-Th2 paradigm" that 
has guided immunologists since the late 1980s states that 
the type of immune response depends on which of the two 
populations prevails in the concurrence of Thl and Th2 
helper cells. For the case of allergy this entails that there 
are populations of allergen-specific T helper cells in both 
allergic and healthy individuals. Yet, in the latter an al- 
lergic resjDonse_is preve nted by the predominance of Thl 
cells (|Romagnanilll997l ). The "hy giene hypothesis" , that 
claims that a hygienic childhood environment increases the 
risk of allergic diseases, can be explained in this framework 



as follows: Due to the reduced exposure to bacterial and 
viral antigens the Thl cells are only insufficiently stimu- 
lated and therefore cannot prevent th e Th2 cells from dom- 



inatin g after exposure to an allergen (lYazdanbakhsh et al 
20021) . 
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In recent years doubts have been raised about the per- 
suasiveness of this explanation. Some studies show, for 
instance, that populations with high rates of helminth in- 
fections are equally protected from allergic diseases, even 
though these infections induce strong Th2-mediated im- 
mune reactions. On the other hand, a considerable in- 
crease in the frequency of type 1 diabetes and other au- 
toimmune diseases, which turn out to be mediated by Thl 
cells, has been observed ( Wills-Karp et al. . 200lh . Thus, 
there seems to be yet another mechanism of regulation 
which is able to prevent the development of unwanted 
immune responses in healthy individuals and whose mal- 
function can lead to either allergic or autoimmune dis- 
ease. A modified version of the hygiene hypothesis known 
as "counter-regulation hypothesis" has been suggested ac- 
cording to which all kinds of infections can possibly pre- 
vent the development of allergic disorders by inducing the 



proliferation of r egula tory T cells (Treg) (jMurphv et al 
2007tlSakagu"chll2000h . 



Different types of regulatory T cells have been iden- 
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tified. The type of Treg cells which seems to be impor- 
tant in the context of allergic diseases is the so-called in- 
duced regulatory T cell (Trl). Cells of that kind produce 
cytokines such as IL10 and TGF-/3 which can suppress 
both Thl and Th2 mediated immune responses and they 
differentiate from naive T cells just as th e other subsets 
( Battaelia etail 120061: rTavlor et all 12006ft . 

Allergen-specific immunotherapy (also known as desen- 
sitization therapy) consists of repeated injections of aller- 
gen or allergen peptides and aims at inducing a state of tol- 
erance in the allergic individual. Even though specific im- 
munotherapy has been carried out for more than one cen- 
tury now, the underlying mechanisms remain poorly un- 
derstood. Within the framework of the Thl-Th2 paradigm 
immunologists assumed that in the course of immunother- 
apy the Th2 mediat ed reaction is "switc hed" to a Thl 
dominated response ([Murphy et all 120071) . More recent 



cells (N) , and the concentration of allergen (A) presented 
by antigen-presenting cells. Figure Q] shows a simplified 
scheme of the T cell interactions that are incorporated into 
the model. Following an injection, allergen is taken up by 



studies indicate, however, that the therapeutic effect is 
mainly caused by an incre ase in the population of allergen- 
speci fic regulatory T cells ( Akdis et al. . 1998 ; Akdis and Akdisl , 

The therapy is performed in practice by starting with 
very small, innocuous injections which are subsequently 
increased until a maximum dose is reached. After that, 
during the maintenance phase, this dose is administered 
once every four weeks over a period of 3-5 years. There 
exist different protocols for the initial part of the treatment 
which differ in the period of time in which they reach the 
maximum dose. In conventional therapies it takes about 
two months, while in so called "rush proto cols" the mainte 
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nanc e dose is reached after only one week (jBousauet et al 
1998ft . 



If we assume that the therapeutic effect of immunother- 
apy is mainly due to a change of T cell equilibrium which 
involves only a small number of cell types, it should be 
possible to capture it within a mathematical model. A 
model using nonlinear differential equations describing the 
dynamics o f Thl- Th2-Interactions has been intr oduced in 
Behn et al.l (2001) and was further investigated in lRichter et al 
(|2002ft and lvoeel and Behnl (|2007t ). On the following pages 
we will present an extended version of the model that takes 
into account the influence of the population of allergen- 
specific regulatory T cells. After motivating the set of 
equations that defines our model and explaining the oc- 
curring parameters, we will investigate the simple case of 
periodic injections. By making use of the stroboscopic 
map, we will already be able to anticipate the qualitative 
features of realistic therapies. Finally, we will show that 
we can simulate different therapy protocols provided that 
the initial conditions are chosen in the right way. 

2. The Model 

Our model consists of a set of nonlinear differential 
equations describing the temporal behavior of five vari- 
ables: the concentrations of Thl, Th2 and Trl cells (Ti, 
Tz, T r respectively), the concentration of naive T helper 



Figure 1: Simplified scheme of T cell interaction in response to al- 
lergen encounter. Allergen presented by an antigen-presenting cells 
(APC) activates naive T helper cells (JV) which leads to their subse- 
quent differentiation to either Thl, Th2, or Treg cells (Ti, T2, and 
T r respectively). Thl and Th2 cells suppress each other and support 
their own proliferation respectively via their cytokines (IF, IL). Treg 
cells suppress Thl and Th2 cells. Dominance of Th2 cells leads to 
allergic reaction. 

an antigen-prcscnting cell (APC) and presented to naive 
T helper cells. Upon activation, these naive cells can dif- 
ferentiate into Thl, Th2 or Treg cells. Via their cytokines 
(IF, IL respectively) activated cells can exert autocrine 
action on their own population and suppress proliferation 
of the other. Thl and Th2 cells suppress each other re- 
spectively, whereas Treg cells suppress both Thl and Th2 
cells while they themselves are not suppressed. The asym- 
metric way in which the populations of Thl and Th2 cells 

I' nteract is adopted from the previous version of the model 
Behn et all 12001 ). Our new attempt to describe the con- 
currence of T cells leads to the following set of equations: 
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We shall now explain the form of the equations as well 
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as the occurring parameters. Looking at equations ([T])- 
we find that the specific T cell populations only grow 
to substantial sizes if allergen is presented. In the absence 
of such a stimulus most cells die off. All T cells (including 
the naive cells) are assumed to have the same half life and 
consequently all populations decay at the same rate. The 
system is already rescaled to dimensionless units, in partic- 
ular the time is measured in units of the half life of T cells. 
Naive cells are produced at a constant rate a, whereas the 
generation of Thl, Th2, and Treg cells is proportional to 
the concentration of naive cells, the concentration of pre- 
sented allergen, as well as to the concentration of their 
respective cytokines (autocrine stimulation). As the cy- 
tokines are degraded fast compared to the half life of cells, 
the concentration of cytokines produced by a T cell sub- 
population can be regarded as proportional to the size of 
that population itself. For that reason, the cytokines do 
not explicitly appear in the equations. The parameter c 
accounts for a small background of cytokines arising from 
other processes of the immune system. It is assumed to be 
equal for the three subsets of differentiated helper cells and 
to be constant over time. Its mathematical role consists 
in initially driving the system away from the trivial state, 
where all T cell concentrations are zero. Suppression is 
modeled by factors of the form 1/(1 + x) where x stands 
for the concentration of cytokines produced by the sup- 
pressing population. In these factors the small cytokine 
background is neglected. Finally, equation §5§ states that 
the presented allergen is degraded proportionally to the 
total concentration of specific T cells. 

The parameter v determines how many differentiated 
T cells arise from one naive cell, <fi and x account for dif- 
ferences in the autocrine action of the three subsets. The 
strength of suppression is regulated by the parameters /ii, 
fi2, and /x r respectively. 

A more detailed derivation for the case of the Thl-Th2 



In the same way, for the case of TijT r we find 



model can be found in lRichter et al.1 (|2002l ). Furthermore, 
it is argued there that the parameters </>, and /12, which 
already occur in the old version of the model, have to sat- 
isfy the conditions </> > 1 and [i\ > fJ-2- In the next step 
we will try to find analogous conditions for the choice of 
X and fi r . To this end, we will turn our attention to the 
development in time of the ratios T±/T r and TijT r . To 
keep the following calculations simple, we use the approxi- 
mation c ~ 0. It can be shown that the conclusions drawn 
also hold for the case of small but nonvanishing c. 
It follows from © and (gj) that 
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Setting this expression equal to zero yields 
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Equations © and (J5J can only be satisfied for positive cell 
concentrations if we set x < 1 an< i X < 4>- Otherwise the 
Treg cells will always dominate over the other two subsets 
which makes it impossible to simulate any allergic reaction 
at all. Provided that x < 1 < 0, we find that above a 
threshold given by 
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(9) 



the Treg cells have a higher growth rate than the two other 
populations. This threshold is independent of the concen- 
trations of Thl and Th2 cells. If its value is set too high, 
the Tregs will never be able to compete and there will be 
no successful therapy. A very low threshold on the other 
hand will make them too dominant. Therefore, the rela- 
tion given by ([9]) can lead us to a reasonable choice of fi r . 
For numerical si mulations, we will alw ays use x = 0.8 and 
fj, r — 0.25. From Richter et al. ( 20021 ) we adopt the choice 
of the remaining parameters: a = 10, v = 8, x — 1-02, 
fix = 0.2, H2 = 0.1, and c = 10~ 4 . 

Equations ((U)-© constitute an autonomous dynami- 
cal system, but this only holds because we have not yet 
considered how the allergen is taken up by the organism. 
In immunotherapy the allergen enters the body via subcu- 
taneous injections. In our model an injection of allergen 
at a given time t is modeled by changing the allergen con- 
centration instantaneously from A(t) to A{t) + D, where 
D > specifies the dose administered. After an injec- 
tion the three T cell populations expand by several orders 
of magnitude. However, after a short time (compared to 
the half life of T cells) the allergen has been degraded 
completely (A ps 0) and the populations will not grow any 
longer. It follows from equations ©-(HI that subsequently 
their concentrations will drop exponentially As the half 
life is the same for Thl, Th2, and Treg cells, the ratios of 
concentrations will by then have rea ched a constant value. 

According to lAkdis et al.1 (|2004j ). it is in particular the 
balance of Th2 and Treg cells that is decisive as to whether 
there will be an allergic reaction or not. In our model we 
can directly compare the initial value of the ratio T2/T,, 
to the constant value that is reached after an allergen en- 
counter. We will therefore call an immune response to a 
given dose D allergic if this ratio has increased compared 
to its initial value, or, mathematically speaking, if 
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Starting from the assumption that the naive cells are in 
their stationary state (N = a) before the allergen en- 
counter and that there is no allergen left from previous 
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encounters (A = 0), the type of reaction only depends on 
the initial conditions T°, T°, and on the allergen dose 
D. 

Of course, this description of allergen administration is 
highly idealized and it only approximates the case when 
allergen is taken up in an injection-like fashion, which ap- 
plies for example to insect stings. Allergic reactions to 
pollen or house dust mite, however, are more complicated 
as the allergen is taken up continuously over time. 

3. Fixed Points and Stable Manifolds 

Our next step will be to describe specific immunother- 
apy, that is, administration of repeated injections. The 
mathematically simplest case is that of periodic injections, 
which means giving the same dose D repeatedly at times 
to + n ■ t, n = 0, 1, 2,3,.. .. The case of periodic injections 
corresponds to the maintenance phase of allergen specific 
immunotherapy. To investigate it in more detail we will 
us e the stroboscopic map. This concept has been proposed 
m IVoeel and Behnl (120071 ) and it has proven to be an im- 
portant tool because it considerably reduces the complex- 
ity of the system. 

We denote by 9(T° , A°;t) the solution of the above 
system at time t for initial conditions 



(T°, A ) = (T°, T 2 °, T! r °, JV°, A ) , 



(11) 



furthermore, let 0t(T°, A ; t) be the projection of this so- 
lution on the three-dimensional subspace of T cell concen- 
trations. Elements of this space are vectors of the form 
T = {T u T 2 ,T r ) T . The stroboscopic map S T ,r>(T) for pe- 
riod r and allergen dose D is then defined as 



S T ,u(T) =6 t (T,A d ;t) 



(12) 



where Ad = (a,D). Equation (fT2")) maps a vector T of 
cell concentrations on the state the system will be in at 
time t = to + t if an injection of dose D is given at time 
t = t$ where T is taken as initial condition. If r is not too 
small, then again at time to + r we will have N ~ a and 
A m 0. This means that applying the stroboscopic map 
repeatedly will be a good approximation for describing a 
periodic therapy. 

We will now investigate the long term behavior of the 
system when such a periodic therapy is applied. It follows 
from what has been said above that we can attack this 
question by simply looking at repeated applications of the 
stroboscopic map. If only one single injection is given, the 
system will eventually reach the trivial state (a, 0, 0, 0, 0). 
However, the perturbations exerted by the repeated injec- 
tions can create periodic orbits in which the concentration 
of each cell subpopulation keeps oscillating in a uniform 
fashion always reaching the same peak value. These or- 
bits correspond to fixed points of the stroboscopic map. 
Numerical simulations show that for a given period r the 
stroboscopic map has up to three stable fixed points and 



several unstable fixed points. If we look at stroboscopic 
maps for different periods, we find that their fixed points 
lie on continuous lines. In figure [5] these are displayed 
as branches of stable and unstable fixed points. Bifurca- 
tions occur at certain critical periods which means that 
the number of fixed points can change if r is changed. 

Applying the stroboscopic map S Tj d repeatedly will 
drive the system to one of its stable fixed points. In gen- 
eral we find three such fixed points, consequently there 
arc three possible outcomes for the corresponding periodic 
therapy In each of the stable fixed points the concen- 
tration of one of the T cell subsets peaks high while the 
concentrations of the two others remain much lower. Ad- 
ministration of periodic injections will therefore always re- 
sult in one cell type eventually dominating the two others. 
Which one of the three subsets will in the end be success- 
ful crucially depends on the initial state given by a vector 
(T^T°,T°) £ R 3 + . We can therefore subdivide the space 
of T cell concentrations into three regions, each being the 
set of all initial vectors leading to the same therapeutic 
result. These regions are just the domains of attraction 
of the stable fixed points of the stroboscopic map. The 
boundaries between the domains of attraction are consti- 
tuted by the stable manifolds of the unstable fixed points. 
In addition to showing the branches of fixed points for 
varying period r, figure [5] displays these stable manifolds 
for the specific example of t = 4. 




100 



Figure 2: Fixed points and stable manifolds. The set of fixed points 
of the stroboscopic map 1X21 1 . for D = 1 and varying period r, which 
is made up of three different branches is displayed. As an exam- 
ple the stable (•••) and unstable (o) fixed points corresponding to 
t = 4 are shown along with the stable manifolds of the unstable 
fixed points (brown surfaces). T cell concentrations are represented 
logarithmically. 

A different way of subdividing the state space is pro- 
vided by the definition of allergic states given in (fl"0|) . The 
boundary of the set of allergic states (corresponding to a 
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reference dose) is given by 



4. Successful Therapy 



lim 



T 2 (t) = T° 

T r (t) T° 



(13) 



that is, by all the states starting from which a single in- 
jection does not change the long term ratio T2/T,.. This 
boundary also forms a two-dimensional manifold in the 
state space and we will refer to it as the separatrix. 

The goal of specific immunotherapy is to drive the sys- 
tem from an initially allergic state to a tolerant state char- 
acterized by increased generation of regulatory T cells. 
In our model this means approaching the Treg-dominated 
stable fixed point of the stroboscopic map. Thus, the ini- 
tial state must lie in the domain of attraction of this fixed 
point. The crucial question therefore is: Are there allergic 
states in the domain of attraction of the "healthy" fixed 
point? In numerical simulations we can show that this 
is actually the case, provided that the period r is not to 
long. In figure [3] the relevant part of the separatrix is 
shown along with one of the stable manifolds correspond- 
ing to t = 1. The states that lie below the separatrix (i.e. 
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Figure 3: Treatable allergic states. Separatrix (light blue) for D = 1 
and one of the stable manifolds (brown) corresponding to Sjy T with 
t — 1. The area between the two surfaces represents the set of 
allergic states starting from which a successful therapy is possible. 

in the allergic region) but above the stable manifold (i.e. 
in the domain of attraction of the healthy fixed point) are 
the ones that allow for a successful therapy. It turns out 
that this set of treatable allergic states increases if we re- 
duce the period between injections. Also it increases if we 
choose a higher allergen dose. We can explain this by re- 
calling our investigation of the ratios of cell concentrations. 
We had found in ((9|), that above a certain threshold the 
Treg cells will have the highest growth rate independently 
of the concentrations of the other two cell types. There- 
fore, even if in comparison there are less Treg cells, they 
may be able to catch up. Short intervals and high doses 
result in high concentrations of all T cell subsets and from 
this especially the Treg cells will profit. 



So far we have only looked at the simplified case of pe- 
riodic therapies, but in computer simulations we can also 
test protocols as they are used in practice. In these pro- 
tocols both intervals between injections and administered 
doses vary. Thus, we cannot directly apply the results that 
we have achieved in the previous section. But it is never- 
theless still true that high doses and short intervals are in 
favor of the Treg population. The different proto cols that 
are used in medical practice (jRueff et al. , 2000) have in 
common that right before the maintenance phase sets in, 
the maximum dose is given several times in short intervals. 
From our point of view this is the decisive step in therapy. 
It ensures that the system reaches the domain of attraction 
of the stable fixed point that the therapy is aiming at. At 
this stage the Treg cells start overriding the Th2 cells. In 
figure 0] simulations according to a conventional protocol 
and according to a rush-protocol are shown. In the conven- 
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Figure 4: Successful therapy. Development of T cell concentra- 
tions according to a conventional protocol (top) and rush-protocol 
(bottom). Initial concentrations in both cases are given by T = 
(0.002, 0.01, 0.003) T . In the conventional therapy the Treg cells start 
dominating at about t = 20, in the rush therapy already at t = 5. 
This corresponds to the time when the highest dose is administered 
for the fi rst time. The ther apies are simulated according to protocols 
found in iRueff et al, 1 1 1200011 . We assume that the maintenance dose 
of allergen corresponds to D = 1 and that one unit of time roughly 
corresponds to one week. 

tional protocol the maintenance phase sets in after about 
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two months, in the rush protocol after only a few days. 
Consequently the therapeutic effect sets in much earlier in 
rush-protocols, which is in accor dance to what has been 
observed in practice ( Coxl l2008f) . In both protocols the 
same periodic orbit is finally reached because the mainte- 
nance phase is the same for both therapies. Therefore, the 
final result of therapy is the same and independent of the 
protocol. 

In our simulations we also find that the Thl /Th2 ratio 
increases during therapy. At the beginning T2 is clearly 
higher than T±, whereas at the end both cell concentra- 
tions have the same order of magnitude. This might be 
an explanation for the "Th2-Thl switch" that has been 
observed during specific immunotherapy. 



5. Conclusion 

Our model is able to describe allergic reactions and the 
course and outcome of allergen-specific immunotherapy on 
the T cell level. Apart from this work and the models that 
it is directly based on, we find other attempts in the lit- 
erature to explain specific immunotherapy by means of 
mathematical descr i ptions of T cell dynamics, for example 
Fishman and Seed |199Q). To our knowledge, however, 
the present paper is the first one to include regulatory T 
cells. We have shown that the basic mechanisms in aller- 
gic reactions can be explained as a competition between 
Th2 and Treg cells. In the model, Treg responses are fa- 
vored by high allergen doses administered in short time 
intervals. Therefore, the decisive event in immunotherapy 
is the beginning of the maintenance phase. Protocols in 
practice could be improved by optimizing this step. As 
immunologists have not yet fully understood all the reg- 
ulatory mechanisms playing a role in allergic diseases, we 
are sure that our model will need further refinement. Nev- 
ertheless, our investigations already provide general tools 
to model immune reactions with interacting lymphocytes. 
Once the immunological picture of T cell regulation is more 
complete, we will be able to give a more adequate descrip- 
tion of the real system. In particular, it would be inter- 
esting to extend the model further by including recently 
identified T cell subsets, such as Th l7 cells, which are 



found to play a role in allerg ic asthma (jOboki et all 12008 
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